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IV 


Objective 


The  objective  of  this  project  is  to  mathematically  investigate  characteristics  of  geodesics, 
describing  possible  minimal  complexity  paths  in  the  group  manifold  representing  the 
unitary  evolution  associated  with  a  quantum  computation. 


2.  Approach 


In  the  Riemannian  geometry  of  quantum  computation  [1-11),  a  Riemannian  metric  can  be 
chosen  on  the  manifold  of  the  (4^  —  l)-dimensional  Lie  Group  SU{2^)  (special  unitary 
group)  of  n-qubit  unitary  operators  with  unit  determinant  {1-27).  The  traceless 
Hamiltonian  of  a  quantum  computational  system  is  a  tangent  vector  to  a  point  on  the 
group  manifold  of  the  n-qubit  unitary  transformation  describing  the  time  evolution  of  the 
system.  The  Hamiltonian  H  is  an  element  of  the  Lie  algebra  su{2^)  of  traceless  2”  x  2^ 
Hermitian  matrices  {25-27)  and  is  tangent  at  the  n-qubit  unitary  operator  U  to  the 
evolutionary  curve  at  t  =  0.  (Here  and  throughout,  units  are  chosen  such  that 

Planck’s  constant  divided  by  27r  is  h  =  1.) 

The  Riemannian  metric  (inner  product)  (., .)  is  a  positive  dehnite  bilinear  form  (iL,  J) 
dehned  on  tangent  vectors  (Hamiltonians)  H  and  J .  The  n-qubit  Hamiltonian  H  can  be 
divided  into  two  parts  P{H)  and  Q{H),  where  P{H)  contains  only  one  and  two-body 
terms,  and  Q{H)  contains  more  than  two-body  terms  (i).  Thus, 

H  =  P{H)  +  Q{H),  (1) 

in  which  P  and  Q  are  superoperators  acting  on  iP,  and  obey  the  following  relations: 

P  +  Q  =  7  PQ  =  QP  =  0,  P^  =  P,  Q^  =  Q,  (2) 

where  I  is  the  identity. 
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(5) 


where  f  denotes  the  adjoint,  and,  except  for  ao,  they  are  traceless. 


Their  products  are  given  by 
Also, 


Tro-j  =  0,  i  7^  0. 

=  I- 

~  i^ijk(^ki  i)])  ^  7^  0) 


expressed  in  terms  of  the  totally  antisymmetric  Levi-Civita  symbol  with  £123 
using  the  Einstein  sum  convention. 


(6) 

(7) 

1,  and 


An  example  of  equation  1,  in  the  case  of  a  3-qubit  Hamiltonian,  is 
P{H)  =  xiai  ®  I  ®  I  +X2(J2  ®  I  ®  I  +  xsas  ®  I  ®  I 

+  X4I  ®ai®  I  +  X5I  ®  a2®  I  +  xqI  ®  (Ts®  I  +  X7I  ®  I  ®  (7i 
+xsl  ®  I  ®  (J2  +  xgl  ®  I  ®  a3+  xioai  ®a2®  I  +  xnai  ®  I  ®  a2 
+X12I  ®  <Ji®  a2+  Xisa2  ®ai®  I  +  XMa2  ®  I  ®  ai  +  xi^I  ®a2®(yi 

+  xiQai  ®a‘i®I  +  xnai  ®  /  0  0-3  -t-  xi^I  ®  ai  0  <73  -h  xiga^  ®  ai  ®  I 

+X20O'3  ®  I  ®  ai+  X21I  ®(Js®Cri+  X22<X2  ®  0-3®  I  +  X23(T2  ®  I  ® 
+X2aI  ®  0-2®  a3+  X250'3  ®  a2  ®  I  +  X260'3  ®  I  ®  0-2  +  X27I  ®  (T3  ®  a2 

+  X23Ckl  ®(Ti®I-\-  X290'2  ®  (72  ®  I  +  XsoOS  ®  <73  ®  I  +  X3iai  ®  I  ®  ai 

-hX32<72  ®  I  ®<72  +  X33a3  ®  I  ®  a3  +  X34I  ®  ai  ®  ai  +  X35I  ®a2®<72 
+X36I  ®  <73®  a3,  (8) 

in  which  (g)  denotes  the  tensor  product,  and 

Q(II)  =  a;37(Ti  (g)  (72  ®  <73  -h  a;38(7i  (g)  0-3  (g)  (72 

+  X3ga2  ®<7i®a3+  a;40(72  (g)  (73  (g)  (7i  (9) 

0:41(73  (g)  (7i  (g)  (72  0:42(73  ®a2®<7i 

+  0:43(71  ®ai®a2  +  0:44(71  ®a2®<7i 
+  0:45(72  (g)  (7i  (g)  (7i  -h  X4eai  (g)  (7i  (g)  (73 
-f  0:47(71  (g)  (73  (g)  (7i  -b  0:48(73  (g)  (7i  (g)  (7i 
+  0:49(72  ®<72®<7i+  0:50(72  (g)  (7i  (g)  (72 
+  0:51(71  ®a2®<72  +  X52<72  ®  (72  ®  (73 
+  0:53(72  (g)  (73  (g)  (72  -b  0:54(73  (g)  (72  (g)  (72 
+  0:55(73  (g)  (73  (g)  (7i  -b  0:56(73  (g)  (7i  (g)  (73 
+  0:57(71  (gi  (73  (g)  (73  -b  0:58(73  (g)  (73  (g)  (72 
+  0:59(73  (g)  (72  (g)  (73  -b  Xeo<72  ®  (73  (g)  (73 
+  Xeiai  (g)  (7i  (g)  (7i  -b  X62<72  ®<72®<72 

+  X63<73  ®  <73  ®a3.  (10) 
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Here,  all  possible  tensor  products  of  one  and  two-qubit  Pauli  matrix  operators  on  three 
qubits  appear  in  P{H),  and  analogously,  all  possible  tensor  products  of  three-qubit 
operators  appear  in  Q{H).  Tensor  products  including  only  the  identity  are  excluded 
because  the  Hamiltonian  is  taken  to  be  traceless.  Each  of  the  terms  in  equations  8  and  9  is 
an  8  X  8  matrix.  The  various  tensor  products  of  Pauli  matrices  such  as  those  appearing  in 
equations  8  and  9  are  referred  to  as  generalized  Pauli  matrices.  In  the  case  of  an  n-qubit 
Hamiltonian,  there  are  4*^  —  1  possible  tensor  products  (corresponding  to  the  dimension  of 
SU{2^)),  and  each  term  is  a  2^x2^  matrix. 

The  right-invariant  {16-18,  26,  27)  Riemannian  metric  for  tangent  vectors  H  and  J  is 
given  by  (i) 

{H,J)  =  ^HlHP(J)  +  gHQ(J)].  (11) 

Here  g  is  a  large  penalty  parameter  that  taxes  more  than  two-body  terms.  The  length  I  of 
an  evolutionary  path  on  the  SU{2^)  manifold  is  given  by  the  integral  over  time  t  from  an 
initial  time  to  a  hnal  time  tj,  namely, 

b 

1  =  (12) 

u 

and  is  a  measure  of  the  cost,  in  terms  of  quantum  circuit  complexity,  of  applying  a  control 
Hamiltonian  H{t)  along  the  path  (i). 

In  order  to  obtain  the  Levi-Civita  connection  on  the  group  manifold,  one  exploits  the  Lie 
algebra  su{2'^)  associated  with  the  group  S'f/(2").  Because  of  the  right-invariance  of  the 
metric,  if  the  connection  is  calculated  at  the  origin,  the  same  expression  applies  everywhere 
on  the  manifold.  Following  reference  1,  consider  the  unitary  transformation 

U  =  (13) 

in  the  neighborhood  of  the  identity  I  C  SU{2^)  with 

X  =  X  ■  a  =  XaCT.  (14) 

cr 

Equation  14  expresses  symbolically  terms  like  those  in  equations  8  and  9  generalized  to  2"^ 
dimensions.  In  equations  12  and  13,  X  is  dehned  in  terms  of  U  using  the  standard  branch 
of  the  logarithm  with  a  cut  along  the  negative  real  axis.  In  equation  13,  for  the  general 
case  of  n  qubits,  x  represents  the  set  of  real  (4*^  —  1)  coefficients  of  the  generalized  Pauli 
matrices  a,  which  represent  all  of  the  n-fold  tensor  products.  Taking  the  trace  of  equation 
13,  it  follows  that  the  factor  x'^  multiplying  a  particular  term  a  is  given  by 

V  =  bTV(A'(r),  (15) 
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The  right-invariant  metric,  equation  10,  can  be  written  as 

=  (16) 

in  which  the  positive  self-adjoint  superoperator  G  is  given  by 

G  =  P  +  qQ.  (17) 

Using  equations  2  and  16,  it  follows  that 

F  =  G-^  =  P +  q-^Q.  (18) 

A  vector  Y  in  the  group  tangent  space  can  be  written  as 

U  =  (19) 

a 

with  so-called  Pauli  coordinates  Here  a,  as  an  index,  is  used  to  refer  to  a  particular 
tensor  product  appearing  in  the  generalized  Pauli  matrix  a.  This  index  notation,  used 
throughout,  is  a  convenient  abbreviation  for  the  actual  numerical  indices  (e.g.,  in  equation 
9,  the  number  57  appearing  in  X57,  the  coefficient  of  ui  (8)  <73  (8)  <73). 

Next  consider  a  curve  passing  through  the  origin  with  tangent  vector  Y  having 
components  y"  =  /dt.  It  can  be  shown  that  the  covariant  derivative  of  a  right-invariant 
vector  field  Z  along  the  curve  in  the  Hamiltonian  representation  is  given  by  {1,  2) 

(VyZ)  =  ^{|y,z]  +  f(iK,G(z)|  +  |z,G(y)])},  (20) 

Because  of  the  right-invariance  of  the  metric,  equation  19  is  true  everywhere  on  the 
manifold. 

The  Riemann  curvature  on  the  group  manifold  affects  the  behavior  of  geodesics  and  can  be 
obtained  as  follows.  In  the  case  of  a  right-invariant  vector  field  Z,  one  has  after  substituting 

Z  =  5^yT,  Y  =  J2v"a  (21) 

T  a 

in  equation  20, 

=  I  ([c^,  r]  +  F([u,  G(r)]  +  [r,  G(a)])).  (22) 

Next  denote  Sq  as  a  set  containing  only  tensor  products  of  the  identity,  and  S'12  as  the  set 
of  terms  in  the  Hamiltonian  containing  only  one  and  two  body  terms,  that  is 

S-Q  (23) 
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and 


U{/  ®  I  ®  ...(Ji  (g)  L.aj  ® 

(24) 

Evidently  then 

W,G(t)]  =  1 

[cT,  r],  reS'i2US'o 
q[(T,T],  T  ^  S12U  So  ’ 

(25) 

and  therefore 

F([<t,G(t)])  =  | 

f  F([(j,r]),  re  512  us'd 
[  qF{[a,T]),  r  ^  5i2  U  5o 

(26) 

Using  equation  18  in 

equation  26,  one  obtains 

f  ([<t,G(t)])=  . 

f  T  e  5i2  U  So 

[  r^5i2U5o  ’ 

(27) 

where 

= 

1  it  [cr,  r]  =  0, 

=  qx  if  r]  oc  A,  and 

(28) 

and  q\  is  dehned  by 

III 

b 

a  e  5o 
<7  e  5i2 
(T  ^  5o  U  5i2 

(29) 

Equation  27  can  also  be  written  as 


F  {[a,  G{t)])  =  [a,  r] .  (30) 

Next  substituting  equation  30  in  equation  22,  and  using  equation  28,  one  obtains 

V<^r  =  r],  (31) 


where 


(32) 


The  Riemann  curvature  tensor  with  the  inner-product  (metric)  equation  16  is  given  by  {29) 


R{W,  X,  Y,  Z)  =  ( Vw  VxU  -  Vx  - 

7i\\v,x]y^  tt) . 

(33) 

After  substituting  the  vector  helds. 

W  =  X  =  V  = 

,  = 

(34) 

(7  a  r  fi 


equation  33  becomes 

Rpar^  =  {'^p'^aT  -  V a'V pT  -  V*[p,<^]r,  fi)  .  (35) 
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Next,  for  three  right-invariant  vector  helds  X,  Y,  and  Z,  one  has 

0  =  Vr  {X,  Z)  =  (A',  VyZ)  +  (VrA,  Z) , 


(36) 


or 

(X,  VyZ)  =  -  ( VyX,  Z) ,  (37) 

and  substituting  equation  34  in  equation  37,  one  then  has 

(a,  Vrlj)  =  -  ( V^a,  /i) .  (38) 

Then  replacing  the  vector  a  in  equation  38  by  the  vector  Vo-r,  equation  31  (see  equation 
7),  one  has 

{VpVaT,  fi)  =  -  {VaT,  V pjj)  ,  (39) 

and  interchanging  indices  p  and  a,  then 


(V^VpT,  ^  -  (VpT,  V^li)  . 


(40) 


Then  substituting  equations  39  and  40  in  equation  35,  and  interchanging  the  hrst  and 
second  terms,  one  obtains 

Rparfi  {V pT,  V ap)  (V q-T,  V p/i)  (V i\^p^tjYi  h}  •  (^1) 


Also  clearly 


V  iyZ  =  iVyZ^ 


SO  equation  41  can  also  be  written  as 


(42) 


^parfi  1  ^ cr/^)  i  ^ pf^'}  ^  * 


(43) 


Next  substituting  equation  31  in  equation  43,  one  obtains  the  following  useful  form  for  the 
Riemann  curvature  tensor  (i): 


Rparp  (i  [p,  t]  ,  i  [(T, /i] ) 

-CayCp^p  {i[a,  T],i[p,  p]) 

-C[p,a],r{i[i[p,(^],r],p).  (44) 

The  geodesic  equation  on  the  SU{2^)  group  manifold  with  the  Riemannian  metric, 
equation  16,  is  obtained  as  follows.  Consider  a  curve  passing  through  the  origin  with 
tangent  vector  Y  having  components  =  dx^ /dt.  The  covariant  derivative  along  the 
curve  in  the  Hamiltonian  representation  is  given  by  {1,  2) 

dZ  i 

(D.Z)  =  {VyZ)  =  —  +  -(\Y,Z]+F  (|K,  G(Z)|  +  |Z,  G(y )])) .  (45) 

(Note  that  the  term  ^  in  equation  45  does  not  appear  in  equation  20  because  there  the 
vector  held  Z  is  taken  to  be  right  invariant,  in  which  case  ^  =  0.).  Equation  45  is  true  on 
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the  entire  manifold  because  of  the  right-invariance  of  the  metric.  Furthermore,  a  geodesic 
in  the  SU{2^)  manifold  is  a  curve  U{t)  with  tangent  vector  H{t)  parallel  transported  along 
the  curve,  namely, 

DtH  =  0.  (46) 

However,  according  to  equation  45  with  Y  =  Z  =  H,  one  has 

D,H='^  +  +  F([HMm  +  \HMH)])),  (47) 

which  when  substituting  equation  46  becomes  (i) 

^  =  (48) 

One  can  rewrite  equation  48  using  the  dual  L  of  if  (i,  2)  and  equation  18, 


L^G{H)=F-\H), 

(49) 

and  then  noting  that 

§  -  1  Y\H))  -  (^)  . 

(50) 

Thus  substituting  equation  48  in  equation  50,  one  obtains 

(51) 

or 

(52) 

and  again  using  equations  50  and  52  it  becomes 

^  =  -ip,L]  =  i[L,ir]. 

(53) 

Furthermore,  again  using  equation  50  in  equation  53,  one  obtains  the  sought  geodesic 
equation  (i): 

^  =  i[L,F(L)].  (54) 

This  equation  is  a  Lax  equation,  a  well-known  nonlinear  differential  matrix  equation,  and  L 
and  iF{L)  are  Lax  pairs  {30-32).  Some  solutions  to  the  geodesic  equation,  equation  54,  are 
given  in  references  1  and  7. 
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3.  Results 


Jacobi  fields  describe  the  divergence  or  convergence  of  neighboring  geodesics  and  are  useful 
in  determining  conjugate  points.  Conjugate  points  are  points  on  a  geodesic  at  which  the 
Jacobi  field  is  vanishing  without  vanishing  in  between  these  points.  It  is  well-known  that 
past  the  first  conjugate  point,  a  geodesic  ceases  to  be  minimizing  {15,  18).  Jacobi  fields  are 
first  to  be  addressed  here  for  a  general  Riemannian  manifold.  Following  this,  Jacobi  fields 
will  be  specialized  to  the  SU{2^)  group  manifold  germane  to  quantum  computation. 


Consider  a  one-parameter  family  of  geodesics  on  a  generic  Riemannian  manifold, 

x^  =  x^{s,t),  (55) 

in  which  the  parameter  s  distinguishes  a  particular  geodesic  in  the  family,  and  t  is  the 
usual  curve  parameter,  which  can  be  taken  to  be  time.  (In  this  section,  Latin  indices  are 
used  in  the  description  of  the  Riemannian  manifold.  Also,  the  in  equation  55  are  not  to 
be  confused  with  the  of  sections  1  and  2.)  The  Riemannian  geodesic  equation  in  a 
coordinate  representation  is  given  by  {18) 

d'^x^  .  ,  dx’^  dx^  ,  , 

in  which  the  Levi-Civita  connection  is  given  by 


2^^  {^){9km,l{s)  +  ()lm,k{s)  9kl,m{s))  ■,  (5 

for  metric  gij{x{s,  t))  =  gij{s).  (Partial  derivatives  are  used  in  equation  56  to  distinguish 
the  s  from  the  t  dependence.)  The  geodesic  equation,  equation  59,  on  the  SU{2^)  group 
manifold  can  be  shown  to  also  follow  from  equation  56  {1,  2). 


Let  x^{Q,t)  be  the  base  geodesic,  and  define  the  lifted  Jacobi  field  along  the  base  geodesic 
by  (i) 

.P{t)  =  ■^^x^{s,t)\s=o,  (58) 

describing  how  the  base  geodesic  changes  as  the  parameter  s  is  varied.  Using  a  Taylor 
series  expansion,  one  has  for  small  As  in  the  neighborhood  of  the  base  geodesic, 

x^{As,  t)  =  x^{0,  t)  +  AsJ^{t)  +  O(As^).  (59) 


Here  x^{As,t)  satisfies  the  geodesic  equation  with  the  metric  gij{As).  Operating  on  the 
geodesic  equation,  equation  56  with  dg  =  ^  and  substituting  equations  58  and  59,  one 
obtains  for  As  — 0, 


(92  ^  As.P{t) 

—  Lim — - - -+ 

dP  As^o  As 


L„„Ap)A!"+a,ri:, 


ox  ox 

w)|s=0  a,  O, 


+riz(0)  ^ 


As.P{t)\  dx^  dx^  d  ^ 

Lim - - -  — — h  — —  —  Lim 

As^o  As  )  dt  dt  dt^s^f) 


dt  dt 
AsJ\t)  I 
As  )  ’ 


in  which  gij{0)  =  Qij  is  the  base  metric  and  r;[;(0)  =  is  the  base  connection.  Equation 
60  then  becomes 


0  = 


f)x^  3t^ 


(9J^  dx''  dx^  dJ^ 
dt  dt  dt  dt 


Taking  account  of  dummy  indices  summed  over,  it  is  clearly  true  that 

,  dx^  dx^  .  ,  dx’^  dx^ 

1  7^-1-  AU  ^  ^  J  i“  -L  ^  ^  .  J  -  U. 


One  also  has 


dt  dt 
j  dx^  dx'' 


r  +  Vi 


kp  mn 

dx"^  dx'^ 


dt  dt  dt  dt 

Also,  using  the  geodesic  equation,  equation  56,  one  has 


JP  =  0. 


kp  g^2  kp  tq  g^  g^  > 

or  renaming  dummy  indices  on  the  right  hand  side,  it  follows  that 


j'^+TLn 


dx^  dx' 


kp  g^2  ^  qk  il  g^  g^ 

Next  adding  equations  61-63  and  65,  one  obtains 


r  =  0. 


0  =  ^  + 


dt^ 


dt  dt 


r.  ,  X  dx'^  dx' 
+'9sr^;(s)|s=o^T-wT-  +  r 


dt  dt 


kl 


{  dJ'^  dx'  dx'^  dJ' 
I  (9f  dt  ^  dt  dt 


J.J  j.q  dx'  dx'  dx'^  dx^ 

■-L  7/^-L  JU  ^  ^  ^  “hi  t.r)i  ^  ^  ^ 


-Vi 


dt  dt 

dxi  Qx^ 


dt  dt 


r  +  v 


kp  mn  g^  g^ 

dx^  dx'^ 


(61) 


(62) 


(63) 


(64) 


(65) 


kp,m  g^  g^ 


TP  I  yJ  •9'^"  TP  I  yJ  r<?  dx'  dx'  , 


kp  g^2 


qk  il  g^  g^ 


(66) 


or  equivalently, 

d^P{t) 

dt^ 


j  UX  UJy  j  ^  UX  L/X  ^ 

^  kLm  r^,  ^  ^  In^  qb  ^ 


dt  dt 


f^rpk 

-pj  pp  J 


'q  ik  g^  g^ 

dx^  dx'  ^ 

n,n— — j'^  +  T: 


dx^  dx' 


kl 


'^P  dt  dt  "  dt  dt 

dJ'^  dx'  dx'^  dJ' 


+ 


dt  dt  dt  dt 


kl\ 


0|s=c 


dt  dt 

dxk  Qx' 


ry^k 

_  pi  ^  ^  jp 
kp  g^2 

dx^n  Qxk 


JP.  (67) 
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Rearranging  terms  in  equation  67,  then 


d^P{t) 

dt^ 


dx^  Qx^ 


dt  dt 


dxk  Qx^ 


f)rY>k  f)rY>m 

j  p  ux  ux 

dt  dt 


jp 


dx^  Qx^ 


dt  dt  dt  dt 

dx'^  dx^ 


■X 


kp,m 


j  U  X 
kp  g^2 


j  dx^  dJ^  j  dx’-  dX 
dx^  dx^ 


-rhr 


-r  -  an 


dxk  gxl 


dt  dt  "  ^^^kiyS)\s-0  g^  g^ 


Recalling  that  the  Levi-Civita  connection  is  symmetric,  one  has 


ri  _  ri 
^  qp  ■*'  ’ 


and  renaming  dummy  indices,  equation  68  becomes 


a  a 
a  2 


f-rJ  ri  I  ri  r?  ri  ■^P^  jk 

ik,l  ^  il,k  ^  Iq^  ik  ^  kp'-  li)  g^  g^ 

j  dx'^  dx^  J  d'^x’^  J  dx’^  dJ’^ 

kp,m^^J  -  kp-^J  -  kl  g^  g^ 

■  dx^  /dP  (9r™  \ 


dxk  gxl 

dt  dt  ’ 


(68) 


(69) 


(70) 


Using  the  well-known  expression  for  the  covariant  derivative  {22,  29),  it  follows  that 

_  d  fDP\  dx^DJP 

DP  ~  di  \  Dt  J  ^  ^P  dt  Dt 

d  /  dJ^  j  dx^  \  j  dx^  DJP 

di  flF  ^  '■plF''  J  ^  '■p  a  nt  ’ 


(71) 


or 


D^P 

DP 


d^P 


gp  ’^P’^  dt 
dx^  ( dJP 


dx^  dx^  jp  ^  d'^x^ 


dt 


kp  g^2 


+ri 


’^P  dt  V  dt 


dx'^ 


^YP 

mn  g^ 


P 


jp  +  r 


j  dx^  dJP 
’^P  dt  dt 


(72) 


Next  the  well-known  Riemann  curvature  tensor  is  given  by  {29) 


f)j  _  pj 
'^ikl  ~  il,k 


pj  I  pj  pP 
ik,l  ^  kp'-  li 


pj  p9 

'-  Iq'-  ik' 


(73) 


Substituting  equations  79  and  73  in  equation  72,  one  obtains  the  so-called  lifted  Jacobi 
equation  {!): 


D^P 

DP 


+  ^iki 


dx^  dx’’ 
dt  dt 


P  +  ar{;(s)|s=o 


dxk  gxl 

dt  dt 


0. 


(74) 
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This  equation  is  useful  for  investigations  of  the  global  behavior  of  geodesics  and  their 
extrapolation  to  values  of  the  parameter  s  characterizing  neighboring  geodesics  on  the 
Riemannian  manifold  (i). 


If  gij  is  independent  of  s,  one  has 


(75) 


the  last  term  of  equation  74  is  then  vanishing,  and  one  obtains  the  standard  Jacobi 
equation  for  the  Jacobi  vector  {18), 


D'^Jl  .  dx^  dx^ 

_ 1_  _ 

Dt^  dt  dt 

Equation  76  is  also  known  as  the  equation  of  geodesic  deviation  {33,  29),  measuring  the 
local  convergence  or  divergence  of  neighboring  geodesics,  and  it  is  useful  in  the 
determination  of  possible  geodesic  conjugate  points  {18,  1).  (Again  it  is  well  to  recall  that 
conjugate  points  are  points  on  a  geodesic  at  which  the  Jacobi  held  is  vanishing  without 
vanishing  in  between  those  points.  It  is  well-known  that  past  the  hrst  conjugate  point,  a 
geodesic  ceases  to  be  minimizing  {15,  18). 


Next  consider  the  factor  in  the  last  term  of  the  lifted  Jacobi  equation,  equation  74, 


^ki  —  '^sr{;(s)|s=o- 


(77) 


Substituting  equation  57  in  equation  77,  one  has 


III 

.2 

i_'^)i_9km,li^'^)  9lm,,ki^'^) 

^  9kl,mi^'^) 

}  . 

)  |s=0 

(78) 

or  equivalently. 

^ii 

^  dg 

3m{g) 

ds  1 

r,nti  + 

s=0  2 

9lm,k  9kl,7n)'> 

(79) 

in  which  one  dehnes 

9  km  =  dsgkm{s)\s=0- 

(80) 

Using  the  well-known  expression  for  the  covariant  derivative  of  a  second  rank  tensor  {29), 
one  has 

9km-, l  9km, I  9kiXml  9mi^kl-  (^^) 

Then  substituting  equation  81  in  equation  79,  one  obtains 


^ki 


dal'^(s)  1  . 

TmfcZ  +  +  9ki^ml  +  9'mJ' 

|s=o  2 


ds 


~^9lm-,k  +  Qli^lnk  +  Qmi^kl 

9kl]m  dkiJ^lm  9liJ^ km)  i 


(82) 
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and  using  equation  69,  then 


^kl  —  2^  idkmil  9lm-,k  9kl-,m) 

^  ds  |,.o 

Next  noting  that 

(s'^g™.)'  =  (^ir  =  0, 

then 

Multiplying  both  sides  of  equation  85  by  T^,  and  using  equation  80,  one  obtains 

9  9mi^  kl  ^  J  ^mkli 

SO  that  equation  83  reduces  to 

^kl  2^  ^9km-,l  9lm-,k  9kl-,m)' 

Finally  then  combining  equations  74,  77,  and  87,  one  obtains 

D'^J^  ,  dj  9x'-  ^  ,  J  J  -^dx^  dx'- 

+  1^9  l£/fcm;Z  +  9ln,-,k  -  9kl,m)  ■ 

Next  dehne  the  vector  held, 

c=  ^  -^9”^{gL,,  +  gL*  -  g«„)-^^. 


(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 


which  is  independent  of  the  Jacobi  held  Jh  Equivalently,  by  symmetry,  equation  89  can 
also  be  written  as 

—  —  .  (90) 

Substituting  equation  81  in  equation  88,  one  obtains  the  second-order  diherential  equation. 


Dt^ 


dxt  g^i 


=  0, 


(91) 


the  so-called  “lifted  Jacobi  equation”  (i).  Nielsen  and  Dowling  used  the  lifted  Jacobi 
equation,  equation  91,  adapted  to  the  SU{2^)  group  manifold,  to  deform  geodesics  by 
varying  the  penalty  parameter  q  (see  section  4).  This  enabled  them  to  dehne  a  so-called 
“geodesic  derivative”  and  to  numerically  deform  a  geodesic  as  the  penalty  parameter  q  is 
varied  without  changing  the  hxed  values  U  =  1  and  U  =  Uf  of  the  initial  and  hnal  unitary 
transformation  corresponding  to  a  quantum  computation  (i). 
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The  generic  lifted  Jacobi  equation,  equation  91,  can  be  solved.  One  first  rewrites  equation 
72  as 


D'^.P 

Dt^ 


d'^P  dx^  dP  ^  dx^  dx^  ^  d'^x^ 

:r"  r  i  ^  J  +  J-  ^ .  o  J 


dt"^ 


dt  dt  dt  dt 


kp  Q^2 


P)^k 

_i_pi  r_rp  7« 

O^-*-  hri  r^  ,  ■*■  mr)  r^  ,  ^  ? 


kp  -  mn 

and  renaming  dummy  indices  in  the  last  term,  then 


D^P 

Dt^ 


d^P 


dt"^ 


2ri 


dx^\  dP 


^P  dt  /  dt 


dx^  dx^  d'^x^  dx^  dx^ 

r-? _ i_  r-? _ i_  r-?  _ r?  _ 

^P’^  dt  dt  ^P  dt^  dt  ^P  dt 


P, 


or  equivalently 


where 


D'^P  d\P  ^PP 
Al- 


Dt^ 


dt"^ 


dt 


\n=l 


Pp  =  ( 2r' 


dx^ 


=  r 


kp  dt  )' 

dx'm  g^k 


P  ’^P’^  dt  dt  ’ 

=  rf 


and 


Next  equation  91  can  written  as 


p  kp  g^2  ’ 

fjrrk  f)rpm 

(3)  Di  =  ri  ^-rq  ^ 

p  -  kq  g^  mp  g^  ■ 


D'^P 

Dt^ 


+  ^^PlP  +  P  =  0, 


where 


(4)  m  W 

dt  dt  ■ 

Next  substituting  equation  99  in  equation  94,  one  obtains 

(9^  P  d  JP 

Ai^  +  BIP  +  C^  =  0, 


dP 


dt 


where 


Next  dehne  the  column  vectors 


BJ  =  A)b^. 

p  p 

n=l 

j  ^  [p] , 


(92) 


(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 

(100) 

(101) 

(102) 

(103) 
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and  the  matrices 


C  =  [C^]  , 

^  ^  [Ai]  . 

■  4 

B  =  M  =  E  '"'^r  ■ 


Equation  101  then  becomes 


(9^  7  (9  7 


Furthermore,  dehning  the  column  vector 


then  equation  107  is  equivalent  to 


J 

=  d.J 


^  ^  0  / 
dt  —B  —A 


(104) 


(105) 

(106) 


(107) 


(108) 


(109) 


The  homogeneous  part  of  equation  109  with  C  =  0  is  equivalent  to  the  Jacobi  equation, 
equation  76,  and  is  given  by 

^  =  MK„,  (110) 


in  which  the  matrix  M  is  given  by 


M  = 


0  I 
-Bit)  -Ait)  ’ 


(111) 


and  the  time  dependence  of  A  and  B  is  indicated  explicitly.  The  solution  to  the  Jacobi 
equation,  equation  110,  is  given  in  terms  of  the  time-ordered  exponential  [34,  25),  namely. 


(112) 


Ko(t)  =  I  +  J2-  *1-  /  *«T(M(«,)...Af(«„))  A'„(0), 


0  0 


where  T  denotes  the  time  ordering  operator  (not  to  be  confused  with  the  transpose  of  a 
matrix,  appearing  below).  Thus,  equation  112  gives  the  Jacobi  field  and  can  be  expressed 
formally  as 


Koit)  =  Texp  /  dt'Mit')  KoiO), 


(113) 


or  dehning  the  operator 


Ei  =  Texp  /  dt’Mit’)  =I  +  J2-  /  /  dtnTiMiti)...Mitn)). 


(114) 


0  0 
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Equation  113  can  also  be  written  as 


A'„(«)  =  E,K„{0). 


(115) 


It  follows  from  equation  114  that 


dEt 

dt 


+  n  /  dti..  /  dtn-iT{M{ti)...M{tn-i)M{t)) 


n=2 


0  0 


M(t)  +  M(t)  ^  j  dti..  J  dtn-iT{M{ti)...M{tn-i)) 

n=2  Q  Q 


=  M{t)  +  j  dtnl^{M{U)...M{tn)) 

\  0  0 

or  equivalently  then  substituting  equation  114,  one  obtains 

dEt 


dt 


=  M{t)Et. 


The  solution  to  the  inhomogeneous  equation,  equation  109  is  given  by 

t 


!  drE;^ 

0 

.  . 

(lie) 


(117) 


(118) 


This  is  the  lifted  Jacobi  held.  To  see  that  equation  118  solves  the  inhomogeneous  equation, 
equation  109,  one  notes  that  using  equations  118  and  117  one  has 


dK{t)  _  dEt_^^^^_dEt 


dt 


dt 


dt 


irE;'  I 


-  EtE^ 


-1 


0 

C{t) 


=  M(t)EtK{0)  -  M(t)Et  /  drE, 


-1 


C{r) 


Next  substituting  equations  115,  118,  and  111  in  equation  119,  then 
dK{t) 


dt 


=  M{t)EtK{0)  +  M{t)K{t)  -  M{t)EtK{t))  - 


0 

Cit) 


0 

Cit) 


0  I 
-Bit)  -Ait) 


Kit) 


0 

Cit) 


(119) 


(120) 


and  thus  equation  109  is,  in  fact,  satished  by  equation  118. 


The  manifold  of  interest  in  the  present  work  is  the  S'17(2”)  group  manifold.  For  this  case, 
consider  a  base  geodesic  with  coordinates  7'^(g,t)  on  the  SU i2^)  group  manifold  with 
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penalty  parameter  q,  and  a  neighboring  geodesic  with  coordinates  7‘^(g  +  A,  t)  with  penalty 
parameter  g  +  A.  To  hrst  order  in  A,  one  has 

-f^{q  + A,t)  =  Y{q,t)  + AJ^{t),  (121) 


in  which  the  Jacobi  held  coordinates  are  dehned  by 


.r{t)  = 


d-f^{q,t) 


(122) 


The  Hamiltonian  for  a  geodesic  'y{t)  with  penalty  parameter  q  is  given  by 

df 

Then  one  has 

dHq  d  d'j 

dq  dq  dt  ’ 

or  equivalently, 

dHq  d  d'y 

dq  dt  dq  ’ 

and  substituting  equation  122  in  equation  125,  then 

dHq  dJ 


in  which  the  Jacobi  held  J  is 


dq  dt  ’ 


J  =  Ha. 


(123) 


(124) 


(125) 


(126) 


(127) 


The  geodesic  equation  for  the  base  geodesic  with  penalty  parameter  q  is  given  by  equation 
53,  namely. 


where  the  dual  L  is  given  by 


=  t[L,H], 


L  =  G{H). 


The  geodesic  equation  for  the  nearby  geodesic  with  penalty  parameter  g  +  A  is 


where 


and  in  accord  with  equation  17, 


L  =  G{H), 


(128) 


(129) 


(130) 


(131) 


G  =  P+{q  +  A)Q  =  G  +  A—  =  G  +  AG', 

dq 


(132) 
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in  which 


(133) 


dq 

Next  letting  U{t)  and  U{t)  denote  the  geodesics  for  penalty  parameter  q  and  g  +  A, 
respectively,  then  for  small  A  one  expects 


U  =  Ue 


and  it  follows  that  to  hrst  order  in  A, 


dt  dt 


But  according  to  the  Schrodinger  equation,  one  has 


(134) 


(135) 


dt 

(136) 

dU 

—  =  -iHU, 

(137) 

in  equation  135,  one  obtains 

-iUA^, 
dt  ’ 

(138) 

or  substituting  equation  134,  then  equation  138  for  small  A  becomes 

-  iUA^,  (139) 

dt 

or  equivalently,  then  to  order  A, 

-iHU  =  -iHU  -  =  -iHU  -  iUA^.  (140) 

dt  dt  ^  ^ 

Next  multiplying  equation  140  on  the  right  by  f/l  and  noting  that  unitarity  requires 

UU^  =  1,  (141) 

then  equation  140  becomes 

H  =  H  +  AU^U^  (142) 

CLL 

to  hrst  order  in  A.  Next  substituting  equation  131  in  equation  130,  one  obtains 


:((?(//))  =  *[(?(//, //|, 


(143) 
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Then  substituting  equations  132  and  142  in  the  left  side  of  equation  143,  the  left  side 
becomes 

4(^))  =  4((g  +  AG')(//  +  AA')),  (144) 

where 

K  =  (145) 

dt  ^  ^ 

Equivalently,  equation  144  to  first  order  in  A  is 

4(^))  =  jG(H)  +  a4(G'(//)  +  G(K)).  (146) 

Next,  using  equations  132,  142,  and  1.45,  the  right  side  of  equation  143  becomes 

=  i[{G  +  ^G'){H  +  ^K),H  +  (147) 

or  equivalently  to  hrst  order  in  A, 

=  i[G{H),H]  + ^{i[G{H),K]+i[G\H),H]+i[G{K),H]).  (148) 

In  terms  of  the  dual,  equation  129,  equation  148  becomes 

f[G(^,;^]  =  i[L,H]  + A{i[L,K]+i[G\H),H]+i[G{K),H]).  (149) 

Next  substituting  equations  146,  129,  and  149  in  equation  143,  one  obtains 

+  A(G'(4//)  +  G(4  A-))  =  *[L,  H]  +  A(i|L,  A|  +  i\G\H),  H]  +  i|G(A),  //]),  (150) 

and  further  substituting  equation  128  in  equation  150,  one  concludes  that 

G'(4a)  +  G(4a'))  =  'ilA,  A]  +  *[G'(A),  H]  +  *1G(A),  H],  (151) 

Furthermore,  multiplying  equation  151  on  the  left  by  G~^,  one  obtains 

G-^G\j^H)  +  j^K  =  G-\i[L,  K]  +  i[G'{H),  H]  +  i[G{K),  H]).  (152) 

But  according  to  equations  128,  129,  and  18, 

j^H  =  tG-^[L,H]  =  tF{[L,H]),  (153) 

so  that,  using  equation  18,  one  has 

G-^G\j^H)  =  FG'ij^H)  =  iFG'F{[L,  H]).  (154) 

Then  substituting  equations  154  and  18  in  equation  152,  one  obtains 

0  =  iFG'F{[L,  H])  +  j^K+  F{i[K,  L]  +  i[H,  G\H)]  +  i[H,  G{K)]),  (155) 
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or 


(156) 


{)  =  j^K+  F{i[K,  L]  +  i[H,  G{K)]  +  G'F{i[L,  H])  +  i[H,  G'{H)]). 

Equation  156  is  the  lifted  Jacobi  equation  for  penalty  parameter  varied  from  g  to  g  +  A 
(i).  It  is  an  inhomogeneous  first  order  differential  equation  in  K. 

If  G'  =  0,  equation  156  reduces  effectively  to  the  conventional  Jacobi  equation,  assuming 
the  form, 

0  =  Ja-  +  F(i[K,  L]  +  i[H,  G[K)]).  (157) 

Equation  157  can  be  rewritten  as  follows  using  equations  2,  17,  and  18: 

j^K  =  -iF{[P{K)  +  Q{K\  P{H)  +  qQ{H)]  +  [P{H)  +  Q{H\  P{K)  +  qQ{K)]).  (158) 

Expanding  the  commutators,  then 

j^K  =  -tF{[P{K\P{H)]  +  q[P{K\Q{H)] 

+  [Q{K\P{H)]  +  q[Q{K\Q{H)] 

+  [P{H\P{K)]  +  q[P{H\Q{K)] 

+  [Q{Hl  P{K)]  +  q[Q{Hl  Q{K)]).  (159) 

The  hrst  and  fifth  terms  cancel,  and  also  the  forth  and  eighth  terms  cancel,  so  one  obtains 

Ja  =  -*(,  -  l)f  (1P(A),Q(A)]  -  1Q(A),P(A)]),  (160) 

or  equivalently, 

Ja  =-i(q-  1)F([P(A),  Q[H)]  +  |P(A),  Q(A)|),  (161) 

or 

jK  =  ^(q-l)F(\Q(H),P(K)]  -  [P(H),Q{K)])^  (162) 

Solving  equation  162  for  K  yields  the  conventional  Jacobi  field. 

The  inhomogeneous  term  in  equation  156  is  given  by 

G  =  F{G'F{i[L,  H])  +  i[H,  G\H)]).  (163) 

Substituting  equations  133,  49,  17,  and  2  in  equation  163,  one  has 

C  =  FQFi[{P  +  qQ){H),P{H)  +  Q{H)]  +  Fi[H,Q{H)] 

=  FQFi  {[P{H\  Q{H)]  +  q[Q{H\  P{H)])  +  Ft  {[P{H)  +  Q{H\  Q{H)]) 

=  FQFi{l  -  q)[P{Hl  Q{H)]  +  Fi[P{Hl  Q{H)].  (164) 
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Equation  164  can  be  rewritten  as  follows: 


C  =  F{P  +  Q)Ft{l-q)[P{H\Q{H)]  +  Fi[P{H\Q{H)] 

-FPFi{l-q)[P{H),Q{H)].  (165) 

But  using  equations  2  and  18  one  has 

PF  =  P{P  +  -Q)  =  =  P,  (166) 

q 

and  using  equations  2,  17,  and  18,  then  equation  165  becomes 

C  =  F\{l-q  +  F-^)t[P{H),Q{H)]-F-^Pt{l-q)[P{H),Q{H)]) 

=  F\1  -q  +  P  +  qQ-P(l- q))t[P{H),  Q{H)] 

=  F^{1  -q  +  q{P  +  QMP{H),  Q{H)],  (167) 


or  using  equation  2,  then 


C  =  F^P(H),Q(H)]. 


(168) 


This  is  a  useful  form  for  the  inhomogeneous  term  in  the  lifted  Jacobi  equation,  equation 
156  (i). 


Next  combining  equations  156,  163,  and  168,  the  lifted  Jacobi  equation  for  varying  penalty 
parameter  q  is  then  given  by 

Ja'  =  *(,  -  l)f  ([Q(A),P(A)]  -  |P(A),Q(A)|)  -  F^P(H),Q(H)].  (169) 

In  terms  of  the  solution  for  K{t),  the  lifted  Jacobi  held  for  varying  penalty  parameter  can 
hrst  be  written  as 

t 

J(t)  =  J(0)  +  j  dt’^.  (170) 

0 

But  according  to  equations  145  and  141,  one  has 

^  =  (171) 

and  substituting  equation  171  in  equation  170,  one  obtains 

t 

J(t)  =  J(0)  +  j  dt'U\t')K{t')U{t').  (172) 

0 

Next  consider  the  case  in  which  the  Hamiltonian  is  constant  along  a  geodesic.  The  geodesic 
equation  52  then  implies 

[G{H),H]  =  Q.  (173) 
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Also,  using  equations  54,  18,  2,  49,  and  133,  one  has 


0  =  —  =  G-^—=zF[L,F{L)]  =  zF[L,P{L)  +  q-^Q{L)] 

=  tF[L,P{L)  +  q-\l-P){L))] 

=  t{l-q-^)F[L,P{L)] 

=  *(1  -  q-^)F[P{H)  +  qQ{H),  P{P{H)  +  qQ{H)] 

=  i{l-q-^)F[P{H)  +  qQ{HlP{H)] 

=  i{q  -  l)F[Q{Hl  P{H)]  =  ^(1  -  q)F[P{Hl  Q{H)] 

=  i{l-q)F[H,Q{H)] 

=  i{l-q)F[H,G\H)]=Q.  (174) 

It  then  follows  from  equations  156,  173,  and  174  that  if  the  Hamiltonian  is  constant,  then 
again  one  obtains  the  conventional  Jacobi  equation,  equation  157.  Thus  if  the  Hamiltonian 
is  constant,  and  J(0)  =  0  and  dJ{0)/dt  =  0,  then  in  accordance  with  equation  206  below, 
J(t)  is  proportional  to  dJ{0)/dt  and  therefore  J(t)  =  0.  In  this  case,  it  then  follows  that 
the  geodesics  for  the  lifted  Jacobi  equation  for  varying  penalty  parameter  remain  the  same 
as  for  the  conventional  Jacobi  equation  and  are  the  same  for  all  values  of  the  penalty 
parameter  q. 

The  so-called  geodesic  derivative  can  be  used  to  determine  geodesics  that  evolve  from  the 
identity  to  a  chosen  unitary  transformation  U  (i).  In  quantum  computation,  one  generally 
wishes  the  quantum  computation  to  evolve  to  some  hnal  unitary  transformataion  which 
solves  a  given  problem.  One  hrst  chooses  a  Hamiltonian  H{0),  which  produces 
U  =  exp{—iH {0)T  at  some  hxed  time  T  along  the  geodesic  for  penalty  parameter  q  =  1. 
The  parameter  q  can  next  be  varied  to  produce  a  corresponding  change  in  the  initial 
Hamiltonian,  and  this  produces  the  so-called  geodesic  derivative  dHq{0)/dq.  Integration 
then  may  produce  a  geodesic  connecting  the  identity  H(0)  =  I  and  the  chosen  unitary 
transformation  U{T)  for  any  penalty  parameter  q. 

To  proceed  then,  the  general  lifted  Jacobi  equation,  equation  156,  for  varied  penalty 
parameter  can  be  solved.  (It  is  convenient  to  solve  equation  156  directly,  instead  of 
equation  169.)  First  substituting  equation  103  in  equation  156,  one  has 

Ja-  =  -i  F{[K,  L]  +  \H,  G{K)])  -  a  (175) 

The  corresponding  homogeneous  equation  is  then 

=  -t  F{[K,,  L]  +  [H,  G{K,)]),  (176) 
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and  it  can  be  solved  if  it  is  first  recast  in  vectorized  form  {36,  37).  For  any  matrix 


Oil  ai2  ...  flln 

021  0-21  ...  02n 

Oml  Om2  ^mn 

one  defines  the  vectorized  form  of  the  matrix  M  by  the  column  vector, 

VeC  M  [uii.-UttiI,  Cli2***^m2****^lri****^'mri]  i 


(177) 


(178) 


with  each  column  of  the  matrix  M  appearing  beneath  the  previous  one,  arranged  in  a 
column  vector.  If  one  has  a  matrix  equation 


C  =  AX  +  XB  (179) 

for  matrices  A,  B,  C,  and  X,  then  it  can  be  shown  that  {36) 

vec  C  =  [{I®  A)  +  {B^  0  /)]  vec  X  (180) 

It  then  follows  that  the  homogeneous  equation  176  can  be  written  in  vectorized  form  as 
follows: 

^(vec  Ks)  =  -iFvec  [KsL  -  LK^  +  HG{Ks)  -  G{Ks)H] 

=  iFvec  [{LKs  +  Ks{-L))  -  {HG{K,)  +  G{K,){-H))] 

=  iF[(/  ®L-L^  ®  /)vec  Ks-{I®H  -H^  ®  /)vec  G(X,)] 

=  iF[{I  ^L-L'^  ®I)  +  {H^  ®I  -I®  H)G]vec  (X«)] 

=  iAvec  Kg,  (181) 


where 

A  =  F[{I®L-L^®I)  +  {H^®I-I®  H)G]  ,  (182) 

and  F  and  G  are  the  vectorized  forms  of  the  superoperators  F  and  G,  respectively  (i).  For 
example,  the  superoperator  F  acting  on  the  matrix  X  can  clearly  be  written  as 

F(A')  =  AWi?,,  (183) 

j 

for  some  matrices  Aj  and  Bj.  But  one  has  {36) 


vec  AjXBj  =  {BJ  ®  Aj)vec  X, 


(184) 


and  using  equation  184  in  equation  183,  then 


vec  F{X)  =  IFbj  ®  Aj)vec  X, 
j 


(185) 
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and  therefore  the  vectorized  form  F  of  the  superoperator  F  is  given  by 


F  =  0  A, 


(186) 


Evidently  the  solution  to  equation  181  is 


vec  Ks(t)  =  T  exp  i  /  A{t')dt'  1  vec  -h"s(0), 


(187) 


where 


vecKsit)  =  Ktvec  Ks{0), 


Kt  =  T  exp  i  /  A{t')dt 


It  follows  from  equation  189  that 


-Kt  =  iA{t)Kt. 


(188) 


(189) 


(190) 


Here  Kt  is  the  propagator  for  the  homogeneous  form  of  equation  175,  namely,  equation  176. 
The  solution  to  equation  175  is  then  given  by 


K(t)  =  unvec  Kfvec  77(0)  —  Ht  drn^  ^vec  C{r)  1  , 


in  which  unvec  unvectorizes  (i),  namely,  for  a  matrix  M, 


(191) 


unvec  (vec  M)  =  M. 


(192) 


To  see  that  equation  191  solves  equation  175,  one  has 


—K{t)  =  unvec  —Ktvec  77(0)  —  {^Kt)  /  dm^^vec  C{r)  —  ^vec  C{t)  1  .  (193) 


Substituting  equation  130  in  equation  193,  then 


K{t)  =  unvec  (  iA{t)Ktvec  77(0)  —  iA{t)Kt  /  drn^  ^vec  C{r)  —  vec  C{t) 


and  substituting  equation  193  in  equation  194,  then 


(194) 


-77(t)  =  unvec  (iH(t)/T;iVec  77(0)  +iH(t)(vec  77 (t)  —  /D^vec  77(0))  —  vec  C(t)),  (195) 
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or 


(196) 


—K{t)  =  unvec  (i74(t)vec  K[t)  —  vec  C{t)) . 

Substituting  equation  182  in  equation  196,  then 

^K{t)  =  unvec  (iF  [(/  0  L  —  L"’"  0  /)  +  ®  I  —  I  ®  H)G]  veciF(t))  —  C{t).  (197) 

But  according  to  equations  176,  181,  and  182, 

unvec  {iF  [(/  0  L  -  O  /)  +  (77^  O  /  -  /  O  i7)G]  vecK{t))  =  -i  F{[K,  L]  +  [77,  G{K)], 

(198) 

and  then  substituting  equation  198  in  equation  197  one  obtains  equation  175 

7a'  =  -i  F([A,  L]  +  Iff,  G(K)\  -  C,  (199) 

as  required. 

Next,  in  order  to  obtain  the  propagator  of  the  standard  (unlifted)  Jacobi  held  J*,  using 
equation  184,  then  equation  171  in  vectorized  form  is 

vec  ®  77'*^)vec  Kg.  (200) 

Substituting  equation  188  in  equation  200,  then 

vec  =  (U^  ®  U'^)KtMec  77s(0).  (201) 

Next  substituting  equation  171  in  equation  201, 

vec  (j^Js)  =  (U^  ®  U^)Ktwec  (7/(0)^  J,(0)f/(0)t),  (202) 

then  for  U (0)  =  7,  one  has 

vec  (^Js)  =  (U^  ®  U^)nt^ec^Js(0)).  (203) 

Unvectorizing  equation  203,  then 

^ Js  =  unvec  [(7/'^  0  U^)Kt(^ec^Js(0))].  (204) 

But  assuming  Js(0)  =  0,  one  has 

t 

Jg(t)  =  j  (205) 

0 
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and  substituting  equation  204  in  equation  205,  then 


Js{t)  =  j  dt'unvec  [(t/"'"  ®  ?7^)fi;i/(vec  ^  Js(0))].  (206) 

0 

Next  defining  the  propagator  jt  that  generates  the  standard  unlifted  Jacobi  field  at  time  T 
by 

(207) 

then  according  to  equation  206,  one  has 


Jt  =  /  dt'unvec  ®  U^)Kt'vec. 


(208) 


It  follows  from  equations  191,  171,  205  and  the  homogeneous  term  having  the  same  form  as 
equation  207  that  at  time  T  the  solution  to  the  lifted  Jacobi  equation  for  varying  penalty 
parameter  is  given  by 

T 

JiT)  ~  —  f  dtU{ty  (  unvec  Kt  (  /  drK~^vec  C{r)  |  |  U{t).  (209) 


If  J(T)  =  0  in  equation  209,  then 

T 


d 

dt 


J (0)  ~  Jt 


-1 


dtU(ty  unvec  Kt  (  /  ^vec  C(r)  U(t) 


(210) 


Next,  according  to  equation  126,  the  Hamiltonian  H„  for  penalty  parameter  q  is  such  that 


(211) 


and  substituting  equation  210  in  equation  211,  one  obtains  the  so-called  geodesic  derivative 

(1) 

d 


■-1 


dq 


Hqi^)  —  jrp 


dtU{t)'^ unvec  Kt  (  /  dvK^^vec  C{r)  1  U{t) 


(212) 


For  g  =  1,  one  has,  according  to  equations  2,  17,  and  18, 

F  =  P  +  -Q  =  P  +  Q  =  G  =  l,  g  =  l, 

q 

and  then  according  to  equation  48  the  Hamiltonian  is  constant, 

dH 


(213) 


dt 


=  0,  g  =  l. 
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(214) 


Then  equation  168  becomes 


C  =  i[P{H),Q{H)l  g  =  l, 
and  then  because  of  equation  214, 

Therefore  in  equation  212  for  g  =  1  one  has  using  equations  68  and  213, 


(215) 


(216) 


Z  =  Kt  \  /  dvK^^vec  C{r)  =  /  dvK^^ivec  [P{H),Q{H)]  ,  (217) 


and  therefore 


dZ 

dt 


^  I  /  drK,  M  +  1 


*vec  [P{H),Q{H)]. 


(218) 


But  for  g  =  1  in  equation  182,  one  has,  according  to  equations  213,  17,  and  49,  F  =  I, 
G  =  1,  and  L  =  H,  and  therefore 


A=[I®H  ®  I +  H'^  ®  I  -I®  H]  =  Q,  g  =  l. 
Then  substituting  equation  219  in  equation  190,  one  obtains 

dnt  _  Q 

dt 

Also,  according  to  equation  189  for  t  =  0,  one  has 

Kt  =  1. 


Then  substituting  equation  220  in  equation  218,  one  obtains 

^=tveclP(H),Q(H)]. 

Also,  according  to  equations  217  and  221, 

Z{0)  =  0. 


Then  combining  equations  222  and  223,  one  obtains  for  g  =  1, 

Z  =  itYec[P{H),Q{H)],  g  =  l, 
or  using  equations  217  and  224,  then 


Kt  I  /  drn^^Yec  C{r)  =  itvec  [P{H),  Q{H)],  g  =  1- 


(219) 

(220) 

(221) 

(222) 

(223) 

(224) 

(225) 
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Substituting  equation  225  in  equation  212,  then  for  g  =  1,  one  has 


d 

dq 


H,(0) 


■  T 

j  dtU{t)Ht[P{H),Q{H)]U{t)  , 

.0 


For  g  7^  1,  one  has,  substituting  equation  168  in  equation  212, 


g  =  1. 


(226) 


d 

dq 


HM 


dtU{t) 


unvec 


drK;\ec  F^P{H),Q{H)] 


Uit) 


g  ^  1. 
(227) 


Equations  226  and  227  give  the  so-called  geodesic  derivative  (i)  which  is  useful  in  the 
numerical  determination  of  optimal  geodesics  representing  a  particular  quantum 
computation. 


4.  Other  Work 


1.  Solutions  to  the  geodesic  equation  were  addressed,  describing  possible  minimum 
complexity  paths  in  the  SU{2^)  group  manifold  and  representing  the  unitary  evolution 
associated  with  a  quantum  computation  (7). 

2.  A  known  obstruction  to  numerically  solving  the  geodesic  equation  is  the  so-called 
Razborov-Rudich  theorem  (i).  It  is  possible  that  the  obstruction,  if  it  occurs,  can  be 
circumvented  by  introducing  and  tweaking  an  additional  parameter  (other  than  g)  in  the 
Hamiltonian.  This  parameter  could  be  one  multiplying  only  two-body  terms. 


5.  Conclusions 


In  this  work  on  the  Riemannian  geometry  of  quantum  computation,  the  Riemann 
curvature  and  geodesic  equation  were  reviewed,  and  the  Jacobi  equation  and  the  lifted 
Jacobi  equation  on  the  manifold  of  the  S'[/(2”)  group  of  n-qubit  unitary  operators  with 
unit  determinant  were  explicitly  derived  using  the  Lie  algebra  su{2^).  The  Riemann 
curvature  is  given  by  equation  44.  The  geodesic  equation  is  given  by  equation  54.  The 
generic  Jacobi  equation  and  its  solution  are  given  by  equations  76  and  112,  respectively. 
The  generic  lifted  Jacobi  equation  is  given  by  equation  91,  and  the  solution  is  given  by 
equation  118.  The  lifted  Jacobi  equation  on  the  SU{2^)  manifold,  for  a  varying  penalty 
parameter,  is  given  by  equations  156  or  169,  respectively,  and  the  solution  is  given  by 
equations  191  and  168.  Also,  the  geodesic  derivative  is  given  by  equation  212.  These 
equations  are  germane  to  investigations  of  the  global  characteristics  of  geodesic  paths  {15, 
18)  and  minimal  complexity  quantum  circuits  {1,  28,  2). 
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7.  Transitions 


1.  This  work  will  supplement  the  small  U.S.  Army  Research  Laboratory  (ARL)/Sensors 
and  Electron  Devices  Directorate  (SEDD)  Mission  Program  on  Quantum  Computing. 

2.  I  initiated  collaborative  research  with  Dr.  John  Myers  and  Prof.  Tai  Tsun  Wu  of 
Harvard  University. 

3.  I  wrote  10  papers  (references  2-11). 

4.  I  presented  a  American  Mathematical  Society  (AMS)  Short  Course  lecture  on  this  work 
at  American  Mathematical  Society  Meeting  in  Washington  DC. 

5.  I  gave  nine  invited  talks  at  international  meetings. 

6.  I  periodically  informed  a  classihed  government  group  of  progress  on  this  research. 
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No.  of 

Copies  Organization 
1  ADMNSTR 

ELEC  DEFNS  TECHL  INFO  CTR 
ATTN  DTIC  OCP 

8725  JOHN  J  KINGMAN  RD  STE  0944 
FT  BELVOIR  VA  22060-6218 

1  CD  US  ARMY  RSRCH  LAB 

ATTN  RDRLCIMG  T  LANDFRIED 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

26HCS  US  ARMY  RSRCH  LAB 
3  CDS  ATTN  IMNE  ALC  HRR 

MAIL  &  RECORDS  MGMT 

ATTN  RDRL  CIM  L  TECHL  LIB 

ATTN  RDRL  CIM  P  TECHL  PUB 

ATTN  RDRL  SES  E  H  BRANDT  (20  HCS) 

ATTN  RDRL  D  J  MILLER 

ATTN  RDRLD  V  WEISS 

ATTN  RDRLD  R  SKAGGS 

ATTN  RDRL  SE  J  PELLEGRINO 

ATTN  RDRL  SE  J  EICKE 

ATTN  RDRL  R  RAO 

ADELPHI  MD  20783-1197 

TOTAL:  31  (I  ELEC,  4  CDS,  26  HCS) 
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